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Recently the OPERA collaboration reported a measurement of a superluminal speed of muon 
neutrinos traveling through the Earth's crust between their production site at CERN and their 
detection site under Gran Sasso, ~ 730 km away. The measurement was based on the assumption 
that the pulse shape of the neutrinos from the decay of parent mesons produced in proton-target 
collisions is the same as that of the incident protons. Here we argue that the effective column 
density of the target along the beam direction decreases with time during the 10.5 fis duration of 
the proton pulse. This is because of the thermal expansion and expulsion of target material along 
the beam by the energy-momentum deposition during the 10.5 /is pulse. The progressive reduction 
in the effective column density during the pulse decreases the neutrino production rate per incident 
proton. It could have advanced the mean production time of the detected neutrinos relative to that 
calculated from the proton pulse-shape, by an amount comparable to the measured neutrino lead 
time 60.7 ± 6.9(stat) ± 7.4(sys) ns. This explanation implies that the planned measurements by 
OPERA of the speed of neutrinos produced in much shorter (a few ns) pulses, should yield a speed 
consistent with the speed of light in free space. 

PACS numbers: 



I. INTRODUCTION 

The OPERA collaboration has reported re- 
cently [l| a measurement of a superluminal speed 
(v/c - 1) = [2.48 ± 0.28(stat) ± 0.30(sys)] x 1CT 5 of 
muon neutrinos traveling through earth between their 
production site at CERN and their detection site under 
Gran Sasso, ~ 730 km away. However, the OPERA 
measurement of the neutrino speed was not based on the 
time difference between the production and detection 
of individual neutrinos but rather on the measured 
distribution of arrival times of neutrinos at Gran Sasso 
with respect to the waveform of the pulsed protons 
entering the CERN Neutrino to Gran Sasso (CNGS) 
system. Although the distance d from the target where 
the parent meson produces the neutrino in the decay 
tunnel is unknown, this introduces negligible difference 
between the time that a proton enters the CNGS system 
and the time that the produced neutrino enters the 
OPERA detector under the Gran Sasso, 



dt d/2cr/^ ~ r m /2 7m <0.1 ns 



(1) 



for pions of lifetime r m = T n ~ 2 x 10 -8 s and Lorentz 
factor 7,„ = 7^ > 100. 

The neutrino speed that was inferred by OPERA was 
based on the assumption that the probability to produce 
the parent meson in the target is the same for all the 
protons in the pulse, and hence, the temporal profile of 
the neutrino pulse is faithfully represented by temporal 
profile of the proton pulse. However, this may not be the 
case because of the large energy-momentum deposition 
in the graphite target along the beam during the 10.5 /is 
proton pulse. This large energy-momentum deposition in 
the target, as we shall argue below, probably reduces the 
meson production rate for late incident protons during a 
pulse due to thermal expansion and expulsion of target 



material from the beam path. This can advance the mean 
production time of the neutrinos relative to the mean 
time of the proton pulse by the neutrino 60.7 ns lead 
time that was inferred by the OPERA collaboration. 



II. THE NEUTRINO PRODUCTION RATE PER 
INCIDENT PROTON 



Let us first demonstrate that a modest decrease in the 
neutrino production rate per incident proton between the 
beginning and the end of the proton pulse can explain the 
OPERA anomaly. The proton pulse shape has a leading 
edge rise-time of about 0.8 fj,s and a trailing edge fall-time 
of about 0.4 /is and approximately a flat top during the 
9.3 /is in between (see Fig. 11 of[l|). For simplicity, let 
us approximate the proton pulse shape by a rectangular 
pulse of duration At = 10.5 /is. Let us also assume that 
due to the change in the effective column density in the 
graphite target along the beam direction, the neutrino 
production rate per incident proton decreases linearly 
with time during the pulse. Consequently, the neutrino 
pulse shape becomes N(t) = N(0) - (AN/ At) t where 
AJV = N(0) — N(At). 

For a small change, AN <^N(0), the mean time of such 
a neutrino pulse is given by, 
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At 
~2~ 



AJV \ 
6AT(0) J 



(2) 



The observed lead time of 60.7 ns in the mean arrival 
time of the neutrinos requires that the effective column 
density encountered by the proton beam decreases by 
w 7% during the proton pulse. 
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III. THE TARGET RESPOND TO THE BEAM 

It would be presumptuous to calculate in detail the 
changes that take place along the beam path in the CNGS 
target during the 10.5 /xs proton pulse. Thus, we shall 
limit ourselves to simple estimates. 

In the OPERA experiment the beam consists of 400 
GeV protons extracted from the CERN Super Proton 
Synchrotron (SPS). The beam cycle is typically 6 s. 
There are two extractions per cycle, each 10.5 fjs, sep- 
arated by 50 ms. Each extraction has 2.4 x 10 13 protons, 
which carry ~ 1.5 x 10 13 erg. The beam is nearly cylin- 
drical with a diameter of 0.5 mm. 

The target is 2 m long sealed container, which contains 
13 graphite rods 100 mm long each. The first two have 
5 mm diameter and the rest have 4 mm diameter. The 
total weight of the graphite is~41 g (for specific density 
of 2.3 g cm -3 ). The specific heat of graphite is 0.17 
calories/g°K, i.e., 7.1xl0 6 erg/g°K. A very large number 
of disks outside the container cool the target by heat 
exchange with the room temperature during the 6 s cycle. 

Most of the prompt energy loss is through the escape of 
energetic particles which are produced in the proton ini- 
tiated showers within the target. However a significant 
fraction of the energy deposited in the target, (~ 10 12 
erg) does not escape during the 10.5 /is pulse and heats 
the target to a very high temperature. The rising tem- 
perature decreases the target density, in particular along 
the beam direction. Well below the Graphite melting 



temperature (~ 3650 °K), the thermal linear expansion 
coefficient of graphite is ~ 7.9 x 10~ 6 (°K)~ yielding 
~ 4% column density decrease at roughly 2500°K. Al- 
though, the thermal expansion coefficient of the target 
at such high temperatures is not well known, a reduction 
of 7% (see section 2) in the effective column density of 
the graphite target along the beam path at the end of 
the proton pulse seems plausible. The thermal radiation 
from the container and the disks cool the target during 
the 6 s cycle back to its much lower initial temperature. 



IV. CONCLUSIONS 

It has been noticed before (e.g. 0), that the assump- 
tion of the same temporal profile of the neutrino pulse as 
that of the proton pulse could be responsible for the neu- 
trino speed anomaly, which was measured by the OPERA 
collaboration. In this short note we have proposed a 
plausible origin for a difference between the proton and 
neutrino pulse shapes and we showed that it could have 
led to the neutrino speed anomaly reported by OPERA. 
In particular, we predict that once the extracted 10.5 /xs 
proton pulses in the CNGS will be replaced by a few ns 
pulses the neutrino speed anomaly will disappear. 
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